Measurement of the Transverse Polarization of Electrons 
Emitted in Free Neutron Decay 
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Both components of the transverse polarization of electrons (ct^ , o"t 2 ) emitted in the /3-decay of 
polarized, free neutrons have been measured. The T-odd, P-odd correlation coefficient quantifying 
(jt 2 , perpendicular to the neutron polarization and electron momentum, was found to be R = 
0.008±0.015±0.005. This value is consistent with time reversal invariance, and significantly improves 
limits on the relative strength of imaginary scalar couplings in the weak interaction. The value 
obtained for the correlation coefficient associated with ar lt N — 0.056±0.011±0.005, agrees with 
the Standard Model expectation, providing an important sensitivity test of the experimental setup. 

PACS numbers: 24.80.+y, 23.40.Bw, 24.70.+S, 11.30.Er 
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Despite the great success of the Standard Model (SM) 
of elementary particles and their interactions, several im- 
portant questions remain open. One of these is the in- 
complete knowledge of physics of CP-violation, or via the 
CPT theorem, time reversal symmetry violation (TRV). 
The SM with the Cabbibo-Kobayashi-Maskawa (CKM) 
mixing scheme accounts for CP violation discovered in 
kaon [l[ and B-meson systems. It fails, however, to 
explain the basic observation of the dominance of mat- 
ter in the present Universe. The discovery of new CP- 
or T- violating phenomena, especially in systems built of 
light quarks with vanishingly small contributions of CKM 
matrix induced mechanisms, would be a major break- 
through. Nuclear beta decay experiments test these sys- 
tems and free neutron decay plays a particular role: due 
to its simplicity it is free of corrections associated with 
the nuclear and atomic structure. Further, final state in- 
teraction effects, which can mimic T violation, are min- 
imal and can, in addition, be calculated with a relative 
precision better than 1% [J|. 

In this Letter we report on the first experiment search- 
ing for the real and imaginary parts of scalar and tensor 
couplings via the measurement of the transverse polar- 
ization of electrons emitted in the decay of free neutrons. 
There exist very few measurements of this observable in 
general and only two in nuclear beta decays. One of 

them, for the 8 Li system Q, provides the most stringent 
limit on tensor coupling constants of the weak interac- 
tion. 

According to [8| , the decay rate distribution as a func- 
tion of electron energy (E) and momentum (p), from 
polarized neutrons is proportional to: 
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which the electron spin is projected and A is beta decay 
asymmetry parameter. N and R are correlation coeffi- 
cients which, for neutron decay with the SM assumptions 
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-1.27, and allowing for a 



small admixture of scalar and tensor couplings Cs, Ct, 
C' s , C' T , can be expressed as: 
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where J is the neutron spin, <x is a unit vector onto 



where S = (C s + C' s )/C v , T = (C T + C' T )/C A and m is 
the electron mass. The R correlation value vanishes to 
the lowest order within the SM. Including final state in- 
teractions it becomes different from zero, Rfsi ~ 0.0006, 
still below the sensitivity of the present experiment. 
A larger measured value would provide a hint for the 
existence of exotic couplings, and a new source of TRV. 

Applying the Mott polarimetry, both transverse com- 
ponents of the electron polarization can be measured 
simultaneously: ut 2 perpendicular to the decay plane 
spanned by neutron spin and electron momentum, rep- 
resented by R, and oji contained in the decay plane and 
associated with N. The SM value of N is finite and well 
within reach of this experiment. Its determination pro- 
vides an important sensitivity test of the experimental 
apparatus. 

The experiment was performed at the FUNSPIN beam 
line at the neutron source SINQ of the Paul Scherrer In- 
stitute, Villigen, Switzerland. The detailed description 
of the design, operation and performance of the Mott 
polarimeter can be found in Q; here a short overview is 
presented. The reported result comprises independent 
analyses of four data collection periods, featuring dif- 
ferent basic conditions like beam polarization, Mott foil 
thickness and acquired statistics (see Table HI)) . 
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FIG. 1: Schematic top view of the experimental setup. A sam- 
ple projection of an electron V-track event is indicated. 

The Mott polarimeter consists of two identical mod- 
ules, arranged symmetrically on both sides of the neutron 
beam (Fig. [I]). The whole structure was mounted inside 
a large volume dipole magnet providing a homogeneous 
vertical holding field of 0.5 mT within the beam fiducial 
volume. An RF-spin flipper (not shown in Fig.[l} was 
used to reverse the orientation of the neutron beam po- 
larization at regular time intervals, typically every 16 s. 
Going outwards from the beam, each module consists of 
a multi-wire proportional chamber (MWPC) for electron 
tracking, a removable Mott scatterer (1-2 /im Pb layer 
evaporated on a 2.5 /im thick Mylar foil, almost trans- 
parent to electrons from neutron decay) and a scintillator 
hodoscope for electron energy measurement. 

A 1 cm thick plastic scintillator allowed for the electron 
energy reconstruction with 33 keV FWHM resolution at 
500 keV. The asymmetry of the light signal collected at 
both ends of the scintillator slab allowed for the deter- 
mination of the vertical hit position with a resolution 
of about 6 cm, while the segmentation (10 cm) of the 
hodoscope in horizontal direction provided a crude es- 
timate of the z-coordinate. Matching the information 
from the precise track reconstruction in the MWPC with 
that from the scintillator hodoscope considerably reduced 
background and random coincidences. 

A 1.3 m long multi-slit collimator defined the beam 
cross section to 4x16 cm 2 at the entrance of the Mott 
polarimeter. In order to minimize neutron scattering 
and capture, the entire beam volume, from the collima- 
tor to the beam dump, was enclosed in a chamber lined 
with 6 LiF polymer and filled with pure helium at atmo- 
spheric pressure. The total flux of the collimated beam 
was typically about 10 10 neutrons/sec. The beam diver- 
gence was 0.8° in horizontal and 1.5° in vertical direc- 
tion. Thorough investigations of the beam polarization 
performed in a dedicated experiment [10( showed its sub- 
stantial dependence on the position in the beam fiducial 
volume. The average beam polarization necessary for 
the evaluation of the N- and i?-correlation coefficients 
has been extracted from the observed decay asymmetry 
using the precisely known beta decay asymmetry param- 
eter A = -0.1173 ± 0.0013 [llj. This approach automat- 



ically accounts for proper integration over the position 
dependent beam density, its polarization and detector ac- 
ceptance. For this purpose single track events (only one 
reconstructed track segment on the hit scintillator side) 
have been recorded using a dedicated prescaled trigger. 
The main event trigger was used to find V-track candi- 
dates: events with two reconstructed segments on one 
side and one segment accompanied by a scintillator hit 
on the opposite side (Fig. [J). 

To extract the beam polarization P the following asym- 
metries were analyzed: 
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where are experimental, background corrected 

counts of single tracks, sorted in 4 bins of the electron 
velocity f3 and 15 bins of the electron emission angle 7 
with respect to the neutron polarization direction. The 
sign in superscripts reflects the beam polarization direc- 
tion. Background counts, i.e. number of electrons not 
originating from neutron decay, were determined by com- 
paring energy spectra of two event classes: (i) the re- 
constructed electron direction crossed the neutron beam 
("from beam") and (ii) the electron origin was outside 
the neutron beam ( "off beam" ) . The procedure relies on 
the assumption that the spectral shape of the background 
is the same for both event classes, while the character- 
istic neutron /3-decay spectrum with end-point energy of 
782 keV is present only in the "from beam" class. For 
this assumption to hold, the "off beam" range has to be 
carefully chosen in both: inclination angle and extrapo- 
lated origin of the tracks on the opposite detector side 
The validity of this method was verified by compar- 
ing background-corrected energy spectra with simulated 
/3-decay spectra in which energy loss and detector resolu- 
tion were taken into account. Such a comparison is shown 
in Fig. (2k- A similar background subtraction procedure 
was applied for the Mott scattering events (Fig. (2b)- In 
the latter case the modification of the (3 spectrum in- 
duced by the energy dependence of the Mott scattering 
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FIG. 2: Background-corrected experimental energy distribu- 
tions (shaded areas) of (a) the single-track and (b) V-track 
events compared with simulations, (c) Background contribu- 
tion (shaded) to vertex i-coordinate distribution of V-track 
events. The arrow indicates the Mott foil position. 
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cross section is clearly visible. Electronic thresholds are 
not included in the simulation — this is why the mea- 
sured and simulated distributions do not match at the 
low energy side. The average neutron polarization values 
for the four data sets arc collected in Table [Til The low 
polarization for the 2004 data set can be traced to a bug 
in the guiding field found post factum and verified in a 
dedicated experiment. 

To extract the N and R correlation coefficients another 
set of asymmetries was considered: 



A (a) = 



i + (a) — n (a) 
i + (a) + n- (a) ' 



(5) 



where n represent background-corrected experimental 
numbers of counts of V-track events, sorted in 12 bins 
of a, defined as the angle between electron scattering 
and neutron decay planes. In the case of V-track events, 
beside the background discussed previously, events for 
which the scattering took place in the surrounding of the 
Mott-target provide an additional source of background. 
Fig-I2fc presents the distribution of the reconstructed ver- 
tex positions in the x-direction for data collected with 
and without Mott-foil. The distribution clearly peaks at 
the foil position. The "foil-out" distribution has been 
scaled appropriately by a factor deduced from the accu- 
mulated neutron beam. 
It can be shown [§| that: 

A (a) - Pf3AT{a) = PS(a) [NQ(a) + R(3H(a)] , (6) 

where the kinematic factors ^(a), G(a) and Ti.(a) repre- 
sent the average values of the quantities J p, J <r and 
JpXff, respectively, S is the effective analyzing power 
of the electron Mott scattering, known in the literature as 
"Sherman function" , and the bar over a letter indicates 
event-by-event averaging. The term P(3AT accounts for 
the /3-decay asymmetry induced nonuniform illumination 
of the Mott foil. Since the f3 and T are known precisely 
from event-by-event averaging, the uncertainty of this 
term is dominated by the error of the average beam po- 
larization P. To evaluate the influence of this term on 
the final result, the fit with the free parameters R and 
N was repeated with P varied by one standard devia- 
tion. The obtained difference enters the budget of the 
systematic errors and is presented in Table |TJ 

TABLE I: Summary of systematic errors for the 2007 data 
set. 



Source 


8N x 10 4 


SR x 10 4 


term Pf3AT 


5 


23 


effective Sherman function S 


29 


8 


guiding field misalignment 


3 


6 


background subtraction 


46 


53 


dead time variations 


8 


0.3 


Total 


55 


59 



Mean values of the effective analyzing powers as a func- 
tion of the electron energy as well as of scattering and 



incidence angles were calculated using the Geant 4 sim- 
ulation framework [l2T |, following guidelines presented in 
This approach allows accounting properly for 
atomic structure, nuclear size effects as well as for the 
effects introduced by multiple scattering in thick foils. 
The accuracy of these calculations has been verified [n| 
by comparison with two experimental data sets: at low 
(120 keV [3) and high (14 MeV 0) electron energies. 

Mapping of the spin holding magnetic field showed 
small nonuniformities in the beam fiducial volume. These 
were corrected for in the analysis. A residual systematic 
effect (see Table UJ was induced by the uncertainty of the 
field measurements. 

The systematic uncertainty is dominated by effects in- 
troduced by the background subtraction procedure, con- 
nected with the choice of the geometrical cuts defining 
event classes "from-beam" and "off-beam". In order to 
estimate this effect, the cuts were varied in a range lim- 
ited solely by the geometry of the apparatus. 

Since the radio-frequency of the spin flipper was a 
source of small noise in the readout electronics, tiny spin 
flipper correlated dead time variations were observed. 
Their influence on the result was corrected for. The resid- 
ual effect is presented in Table UJ 

The asymmetries as defined in Eqs.S] and E] have been 
calculated for events with energy larger than the neutron 
/?-decay end-point energy and for events originating out- 
side the beam fiducial volume and were found to be con- 
sistent with zero within statistical accuracy. This proves 
that the data were not biased by, for instance, a spin 
flipper related false asymmetry. 

A fit of the experimental asymmetries A, corrected for 
the P(3AT term, to the experimental data set of 2007 is 
shown in Fig. (3) The R and N coefficient values extracted 
in this way from all data sets are listed in Table [TTJ 
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FIG. 3: Left panel: experimental asymmetries A corrected 
for the PfiAT term for the 2007 data set as a function of a 
(defined in text). The solid line illustrates the two-parameter 
(TV, R) least-square fit to the data. Indicated errors are of sta- 
tistical nature. Right panel: geometrical factors T(a), G(ot) 
and 7i(a) for the same data set. Dotted lines are to guide the 
eye only. 

From the approximate symmetry of the detector with 
respect to the transformation a — > — a it follows that 
(3, S and the factors f ', TL are all symmetric while Q is 
antisymmetric function of a (see Fig. [3]). This allows the 
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TABLE II: Summary of results obtained in all data collection periods. Statistical and systematic uncertainties follow the 
experimental values, d is the nominal thickness of the Mott foil, V n represents the total, background corrected number of 
Mott-scattered events, Ek is the average kinetic electron energy for those events and Nsm is the SM value of the N coefficient 
calculated at Ek- Its error comes from the experimental uncertainty of the decay asymmetry parameter A [TTj ] . 



Run 


d (jim) 




PxlO 2 


E K (keV) 


iVsA/xlO 3 


iVxl0 3 (Eq.[7ll 


TVxlO 3 (EqJBJ 


R xlO 3 


2003 


1 


19000 


80.3±1.3±1.6 


331±1.0 ±15 


71±1 


110±108±27 


82±97±31 


-89±143±38 


2004 


1 


74000 


44.2±0.4±1.5 


368±0.5 ±15 


68±1 


144±92±15 


70±86±17 


-117±140±26 


2006 


2 


312000 


80.0±1.0±1.5 


365±0.2 ±10 


68±1 


79±32±7 


86±30±8 


-11±42±9 


2007 


2 


1747000 


77.4±0.2±0.7 


370±0.1 ±10 


68±1 


54±12±5 


51±12±6 


12±16±6 


Total 




2152000 








59±11±4 


56±11±5 


8±15±5 



extraction of the N coefficient from the expression Q: 



N 
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The advantage of this method is that the effect associated 
with the term PfjAT is suppressed by a factor of about 
60 as compared to Eq.[B] The good agreement between 
the N values obtained in both ways enhances confidence 
in the extracted N and R coefficient values. 

In Fig.g] the new results have been included in exclu- 
sion plots containing all experimental information avail- 
able from nuclear and neutron beta decays as surveyed in 
Ref. [l7| ■ The upper plots contain the normalized scalar 
and tensor coupling constants S and T (see Eq.HJ, while 
the lower ones correspond to the helicity projection am- 
plitudes in the leptoquark exchange model, as defined in 
Ref. This is the first determination of the N cor- 

relation coefficient. Although the present accuracy does 
not improve the already strong constraints on the real 
part of the couplings (left panels), the obtained result is 
consistent with the existing data and, in addition, adds 
confidence to the validity of the extraction of R. As to 
the imaginary part (right panels), the new experimental 
value of the R coefficient significantly constrains scalar 
couplings beyond the limits from all previous measure- 
ments. The result is consistent with the SM. 
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FIG. 4: Experimental bounds on the scalar vs. tensor nor- 
malized couplings (upper) and leptoquark exchange helicity 
projection amplitudes (lower panels). T he g rey areas repre- 
sent the information as defined in Ref. while the lines 
represent the limits resulting from the present experiment. 
Solid, dashed and dotted lines correspond to 1-, 2- and 3- 
sigma confidence levels, respectively, in analogy to decreasing 
intensity of the grey areas. 
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